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Motivated by the observation of a disordered spin ground state in the S — 3/2 material 
Bi3Mn40i2N03, we study the ground state properties and excitation spectra of the S = 3/2 (and 
for comparison S = 1/2) bilayer Heisenberg model on the honeycomb lattice, with and without 
frustrating further neighbor interactions. We use series expansions around the Neel state to calcu- 
late properties of the magnetically order phase. Furthermore, series expansions in 1/A = Ji/J±, 
where Ji is an in-plane exchange constant and J± is the exchange constant between the layers are 
used to study properties of the spin singlet phase. For the unfrustrated case, our results for the 
phase transitions are in very good agreement with recent Quantum Monte Carlo studies. We also 
obtain the excitation spectra in the disordered phase and study the change in the critical A when 
frustrating exchange interactions are added to the 5 = 3/2 system and find a rapid suppression of 
the ordered phase with frustration. Implications for the material Bi3Mn40i2N03 are discussed. 

PACS numbers: 74.70.-b,75.10. Jm,75.40.Gb,75.30.Ds 



INTRODUCTION 

The honeycomb lattice is a bipartite but low- 
coordination number two-dimensional lattice. Thus co- 
linear antiferromagnetism is unfrustrated on the lattice. 
However, it is more susceptible to disorder due to its 
low coordination number. Recent interest in antiferro- 
magnetism on the honeycomb lattice comes from many 
directions. First, theoretical studies of quantum spin 
models on the lattice find a rich phase diagram, with 
several colinear, spiral and spin-disordered phases, [l- 14 
Second, a recent Quantum Monte Carlo finding that the 
Hubbard model on the honeycomb lattice shows strong 
evidence for a quantum spin-liquid phase sandwiched 
between the semi-metal phase at small U/t and an or- 
dered antiferromagnetic phase at large U /t, has led to 
much follow-up activity. [1 5(4 19| Third, the discovery of 
graphene and the search for correlated topological insula- 
tor phases have also led to interest in strongly correlated 
electron models on the honeycomb lattice, including spin- 
orbit coupling. [2!" 
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In this paper, our primary motivation comes from the 
bismuth manganese oxynitrate material Bi3Mn40i2N03, 
an 5 = 3/2 antiferromagnet, with a Curie- Weiss constant 
of order 250K, which does not show any long range or- 
der down to 0.4i4r.[22| The material consists of honey- 
comb lattices of S* = 3/2 Mn spins which are separated 
by bismuth and nitrate layers. Two such layers are sep- 
arated by bismuth atoms, forming a bilayer, and these 
bilayers are then separated by significantly larger sepa- 
rations. Thus an appropriate model for such a system is 
a spin-3/2 Heisenberg model on the bilayer honeycomb 
lattice, with an in-plane nearest-neighbor exchange Ji 



and an exchange between the bilayers of Jj_.[23| The bi- 
layer exchange J± has been estimated by an electronic 
structure calculation to be between Ji and 2Ji, whereas 
further neighbor in-plane exchanges, which frustrate the 
system are down compared to Ji by an order of mag- 
nitude (less than 0.2Ji). 23] This system has also been 
studied recently by Quantum Monte Carlo (QMC) sim- 
ulations and bond-operator based mean-field and varia- 
tional theory [23, [111, which our results are compared to 
below. 

We use series expansion methods to study the proper- 
ties of the bilayer honeycomb Heisenberg model., 26, 23| 
Ising series expansions are used to calculate the ground 
state energy and antiferromagnetic order parameters in 
the Neel phase. Dimer series expansions are used to cal- 
culate the ground state energy, triplet energy gap and the 
excitation spectra in the disordered singlet phase. The 
phase transition in this model is known to be in the uni- 
versality class of the 3d classical Heisenberg model. 24 1 
We use this in our series analysis to make our analy- 
sis more accurate. For the nearest neighbor model, we 
find the phase transition to be located at Jx/Ji value 
of 1.66 ± 0.01 for S = 1/2 an 9.34 ± .20 for S = 3/2. 
These results are in very good agreement with the values 
reported in the QMC study of 1.645(1) and 9.194(3) for 
the two cases respectively. The bond-operator theorvpsj 
is significantly less accurate. 

Since the material Bi3Mn40i2N03 is unlikely to have 
J±/Ji « 9, frustration in the plane must play a role in 
disordering the system. Thus, we also study the model 
with frustration and discuss its possible relevance to the 
materials. We also present results for the excitation spec- 
tra in the disordered phase, which could be helpful in 
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further quantifying the experimental system. 

SERIES EXPANSIONS FOR THE HONEYCOMB 
BILAYER MODELS 

We consider the Heisenberg model on the bilayer hon- 
eycomb lattice with Hamiltonian 

(1) 

Hear a = 1,2 denotes the spins in the two layers. The 
first sum runs over nearest-neighbors in the honeycomb 
planes and the second sum is over the second neighbors 
in the honeycomb planes. The third sum is over the 
neighboring spins between the two bilayers. 

For the unfrustrated models (J2 = 0), we have calcu- 
lated Ising expansions for the ground state energy and 
sublattice magnetization to 14-th order for S" = 1/2 and 
to 12-th order for S — 3/2. The dimer expansions in 
1/A = Ji/J± are calculated for the ground state energy 
to order 12 for 5" = 1/2 and to order 9 for S* = 3/2. The 
energy gap series is calculated to order 10 for S — 1/2 
and to order 6 for S = 3/2. For J2/J1 non-zero, the 
energy gap series is calculated in powers of 1/A to order 
6 for = 3/2. The series can be made available upon 
request. 

In Fig. 1 (a) and (b), we show plots of the ground 
state energy for J2 = 0. For both 5* = 1/2 and S = 
3/2, the results from Ising and dimer expansions join 
smoothly as expected for a second order phase transition. 
However, the energy is not the best quantity to determine 
the location of the phase transition. 

In Fig. 2 (a) and (b) , we show the Neel order parame- 
ters for the S — 1/2 and = 3/2 cases. We use a square- 
root transformation[28] to remove the singularity before 
analyzing the series by Fade approximants. Our analysis 
is not accurate close to the transition where the magne- 
tization vanishes. The uncertainties clearly become large 
for the shorter S" = 3/2 series for A > 8 as the transi- 
tion is approached. Beyond A = 9, the rapid decrease is 
sketched in the plot by a dashed line. The order param- 
eter is also not the best way to get a precise estimate of 
the transition point, because the series are not directly in 
the variable 1/A, in which a power-law singularity, with 
an exponent /? is expected. Rather, for every A the series 
is analyzed in the anisotropy variable to calculate the or- 
der parameter for the Heisenberg model. It is difficult 
to enforce the correct power-law singularity in such an 
analysis. 

The triplet excitation spectra for S* = 1/2 and 
S — 3/2 along symmetry lines of the Brillouin 
zone (r(0,0),i^(|f,0),M(|f,^),r(0,0)) are shown in 
Fig. 3 and Fig. 4 respectively. In the dimerized phase, 
there are two branches of the triplet spectrum at every k 
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FIG. 1: Ground state energy as a function of A = J±/Ji for 
(a) 5=1/2 and (b) S = 3/2 models. 



corresponding to the two atoms per unit cell in the hon- 
eycomb lattice. They become degenerate at the M point. 
The spectral gap approaches zero as the transition to the 
Neel phase is approached. 

The most accurate way to determine the critical A in 
our study is by analyzing the energy gap series. Here, 
we use the knowledge that this series has a singularity at 
the critical A with a critical exponent v given by the 3d 
Heisenberg universal value v = 0.71 [2^. Thus, we first 
raise the series to a power l/i/ ^ 1.408, and then study 
it by Fade approximants to see where it goes to zero. 
The Fade approximants show excellent internal consis- 
tency for the longer S" = 1/2 series. Three different ap- 
proximants are shown in Fig. 5a and they are almost 
indistinguishable from each other. Thus the uncertainty 
comes primarily from varying the value of v. We esti- 
mate 1/Ac = 0.6014±0.002, which translates to J_l/ Ji = 
1.66 ± 0.02. For S = 3/2, the shorter series though less 
accurate than for s — 1/2 is still very accurate. Two 
different approximants are shown in Fig. 5b, and one can 
see that they begin to deviate a little close to the tran- 
sition. In this case, we estimate 1/Ac = 0.107 ± 0.002, 
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FIG. 6: Phase Diagram for the S = 3/2 model with frustra- 
tion in the honeycomb layers. 



which translates to J±/Ji = 9.3 ± 0.2. Both these re- 
sults are in very good agreement with the QMC study, 
which gives JjJi = 1.645(1) for S = 1/2 and 9.194(3) 
for S = 3/2 respectively. [2J| 



IN-PLANE FRUSTRATION AND THE BISMUTH 
MANGANESE OXYNITRATE MATERIALS 

Since the spin-3/2 material Bi3Mn40i2N03 is unlikely 
to have J± values an order of magnitude larger than 
in plane couplings, it is clear that frustration must be 
present to explain the absence of Neel order. To study 
this we add a frustrating second neighbor interaction J2 
to our study of the S = 3/2 model. We find that the ad- 
dition of frustration rapidly decreases the critical value of 
J±/Ji- The estimated phase boundary is shown in Fig. 6. 
For J2/J1 = 0.1 this value is reduced to 5.5 ± 1, which is 
still significantly larger than the estimated ratio between 
1 and 2 for the material in the electronic structure cal- 
culation. We find that a value of J2/J1 larger than 0.15 
can lead to a disorder with J±/Ji less than 2. These 
numbers are within the range of values estimated for the 
material from electronic strusture calculations. [23| 

At the classical level, a single layer Heisenberg model 
loses Neel order at J2/J1 = 1/6-0 This is indicated 
by X in Fig. 6. However, this ratio is somewhat en- 
hanced in the presence of quantum fluctuations. Since, 
this enhancement in the critical J2/J1 is roughly 0.03 
for the S = 1/2 case, [13] it should be much smaller for 
the S" = 3/2 case. The critical J2/J1 should be further 
enhanced by the addition of a small J± in the bilayer, 
since a small bilayer coupling makes the system more or- 
dered. Thus, the phase diagram in the J2/J1 and J±/Ji 
plane has a slightly renetrant character for small J±/Ji. 
This has the effect that beyond J2/J1 = 0.17, the phase 
boundary becomes close to vertical and our series analy- 



sis is no longer accurate to locate it precisely. We sketch 
this phase boundary in Fig. 5 by dashed lines. We have 
analyzed the series at J2/J1 = 0.2 and it shows no con- 
sistent point where the gap vanishes. 

If the material Bi3Mn40i2N03 has J2/J1 < 0.17, then 
the singlet phase must arise from strong bilayer cou- 
pling and the material is almost certainly in the spin- 
disordered phase that is adiabatically related to the prod- 
uct singlet phase at large J±. In this phase our calcula- 
tion of the excitation spectra should be accurate. Neu- 
tron scattering on the material can help determine more 
precise exchange parameters. 

However, it is possible that the material has J2/J1 > 
0.17 and in this case depending on how large J±/Ji is, 
a magnetically disordered phase may or may not be adi- 
abatically related to the product singlet phase found at 
large J±. At J±/Ji < 2, there maybe a phase transition 
from one singlet phase to another. Whether a new spin- 
liquid phases arises in the spin-3/2 model at small J±/ Ji 
with frustration, and how far down in J±/ Ji the product 
singlet phase continues remains an open question, which 
deserves further study. If there is a frustration domi- 
nated spin-liquid phase for S = 3/2 at small J± values, 
and that is the appropriate phase for the material, that 
would make the material much more interesting. In this 
case, neutron scattering may show an absence of triplon- 
like excitations. Further experimental study is needed 
before more conclusions can be drawn. 

CONCLUSIONS 

In conclusion, in this paper we have used series ex- 
pansion methods to study the ground state properties of 
S" = 1/2 and 5 = 3/2 bilayer honeycomb lattice Heisen- 
berg models. We find that an explanation for the mate- 
rial Bi3Mn40i2N03 requires significant frustration in the 
honeycomb planes. We also present results for the evolu- 
tion of the spectra with A = J±/Ji in the spin-disordered 
phase that is adiabatically related to the product singlet 
phase at large J± . Further study of the experimental sys- 
tem, especially its triplet excitation spectra is necessary 
before further conclusions can be drawn as to whether 
the material is in the product singlet phase or in a novel 
spin-liquid phase. 
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